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ABSTRACT 

Dark energy observations may be explained within general relativity using an inhomogeneous 
Hubble-scale depression in the matter density and accompanying curvature, which evolves 
naturally out of an Einstein-de Sitter (EdS) model. We present a simple parameterization 
of a void which can reproduce concordance model distances to arbitrary accuracy, but can 
parameterize away from this to give a smooth density profile everywhere. We show how the 
Hubble constant is not just a nuisance parameter in inhomogeneous models because it affects 
the shape of the distance-redshift relation. Independent Hubble-rate data from age estimates 
can in principle serve to break the degeneracy between concordance and void models, but 
the data is not yet able to achieve this. Using the latest Constitution supernova dataset we 
show that robust limits can be placed on the size of a void which is roughly independent of 
its shape. However, the sharpness of the profile at the origin cannot be well constrained due 
to supernova being dominated by peculiar velocities in the local universe. We illustrate our 
results using some recently proposed diagnostics for the Friedmann models. 
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1 INTRODUCTION 

An odd explanation for the dark energy problem in cosmology 
is one where the underlying geometry of the universe is sig- 
nificantly inhomogeneous on Hubble scales, and not homoge- 
neous as the standard model assumes. These models are pos- 
sible because we have direct access only to data on our null- 
cone and so can't disentangle temporal evolution in the scale fac- 
tor from radial variations. Such explanations are ungainly com- 
pared with standard cosmology because they revoke the Coper- 
nican principle, placing us at or very near the centre of the uni- 
verse. Perhaps this is just because the models used - Lemaitre- 
Tolm an-Bondi (LTB) or Szekeres to date /Moffat & Tatarski' 
1992. 1994; M ustapha et al. 1997; Sugiura et al. 1999; Alnes et al. 
2006l;ICele rier 2006;'Vanderveld et al. 2006; Alexander et al. 200J; 
CelerieJ l2007i; iClifton et al.i ,2008.; Garcia-Bellido & Haugb0S 



2008al; iGarcia-Bellido & Hatigb0llell2OO8bl; llshaketalj 20081 ; 
Yooetal.1 1200 8'; 'Zibin et al.' '2008'; 'B oleiko & Wvithd |2009| ; 

200^ - are very simplistic 



Celerier et al. 



2009; Ouercellini et al. 
descriptions of inhomogeneity, and more elaborate inhomogeneous 
ones will be able to satisfy some version of the Copemican prin- 
ciple (CP) yet satisfy observational constraints on isotropy (e.g. , 
a Swiss-Cheese model jMarra et al]|2007l ; iBiswas & Notarill2008l ; 



iMarra et al.|[2008l) or something like that)Q We may instead think 
of these models as smoothing all observables over the sky and so 
compressing all inhomogeneities into one or two radial degrees 
of freedom (d.o.f.) centred about us - and so we needn't think 
then as ourselves 'at the centre of the universe' in the standard 
way0 Whatever the interpretation, such models are at the toy stage, 
and have not been developed to any sophistication beyond under- 
standing the background dynamics, and observational relations; 
in particular, perturbation theory and structur e formation is more- 
or-less unexplored, though this is changi ng l lTomitij|T997l ; IZibiiJ 
l2008l ; IClarkson. Clifton & Febmarvll200a) . They should, however, 
be taken seriously because we don't yet have a physical explana- 
tion for dark energy in any other form. Indeed, one could argue that 
these models are in fact the most conservative explanation for dark 
energy, as no new physics needs to be introduced. 

Regardless of the details, these models raise an important 
question for cosmology, particularly so in the light of the dark en- 
ergy problem; can we test the Copemican principle? While many 
would argue that it has effectively been done via the success and ac- 
curacy of the standard model, until we have a physical explanation 
for dark energy (and the inflaton for that matter) we are open to the 
accusation of having only phenomenological descriptions of two 
key observations in cosmology: accelerating expansion and scale- 
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1 See lClarkson & BarretJ 0993) and iBarrett & Clarksoij )2000l) for exam- 
ples of spacetimes which, although unrealistic, are globally inhomogeneous 
yet satisfy a version of the cosmological principle. 
^ As argued by K. Bolejko and M.-N. Celerier, private communication. 
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invariant initial conditions. Only a handful of tests have actually 
been proposed: the Goodman-Caldwell-Stebbins test, which looks 
at the Cosm ic Microwave Background (CMB) inside our past light- 
cone dOoodman 1995 ; Cal dwell & Stebbins 200S); and the 'curva- 
ture test', which checks if today's value of the curvature parameter, 
r2fe, given by 



0.k = 



{H{z)[{l + z)dA{z)l,Y 



(1) 



[Ho{l + z)dA{z)f 

yields the same answer regardless of the r edshift of measure- 
ment jClarkson et"al]|2008l : lUzan et alj |2008|) . as it must in any 
Friedmann-Lemaitre-Robertson- Walker (FLRW) universe. A cru- 
cial issue for this test is that it requires two independent measure- 
ments: one for distances, from, e.g. Type la Supernovae (SNla), 
and one for the Hubble rate, from, e.g. Baryon Acoustic Oscil- 
lations (BAO) or age estimates. However, recognising the diffi- 
culties posed by the requirement of independent observables that 
these tests require, it has recently been proposed in stead that we 
can t est the Copernican principle from SNIa alone jClifton et al.l 
l2008h . thereby making the process much simpler than these other 
methods. If one dema nds that the vo i d is su itably smooth at the cen- 
tre, then, as argued in lClifton et al.l ( |2008| ). this leads to inevitable 
differences from the A cold dark matter (ACDM) distance modu- 
lus which can be detected with future SNIa observations such as 
the SuperNova Acceleration Probe (SNAP). The essence of this ar- 
gument is as follows. A generic LTB model has two radial func- 
tional d.o.f.. One of these is the bang time function, which con- 
trols any inhomogeneity of the big bang surface; the other may 
be chosen to be the radial curvature today. These functions can 
be chosen such that the observer at the centre observes distances 
and H (z) (or number-counts) exactly as one w ould in a ACDM 
model jMustapha et al .ll997l;ICelerier et al.l2009l) . (In this scenario 
the curvature test would have to be performed using both radial and 
tangential Hubble rates - see below for definitions.) However, if the 
bang time fu nction is n ot constant then this excites modes which 
are decaying Jsii^l 19771) ; hence, if this is significant enough a fac- 
tor at late times to affect these observables the models would have 
to be outrageously inhomogeneous at early times - this doesn't 
rule them out a priori, but does mean such models would have 
to have their early time behaviour rigourously examined. One can 
argue on the basis of this that a 'realistic' LTB model has only 
one true d.o.f.. Now, if this d.o.f. is chosen to reproduce the dis- 
tance modulus of a ACDM model, the void profile must be very 
spiky, and have a discontinuity at the origi n. Although there 
is som e deba te as to whether this i s unrealistic jVanderveld et akl 
I2OO6I) or not jKrasinski et alj2009t) . lciifton et alj ( I2OO8I) show that 
this, when combined with an assumption of asymptotic flatness, 
can be critical in deciding if SNIa contain enough information to 
test the Copernican principle: the distance modulus must differ sig- 
nificantly from the ACDM one. 

In this paper we investigate this issue using the latest SNIa 
dataset, which includes SNIa at very low redshift (^min ~ 0.015). 
The void profile can be made sufficiently differentiable by shaving 
off the spike in the void profile at any radius. Below z = Zmin, 
objects are not entirely in the Hubble flow and thus measurements 
are dominated by peculiar velocities, meaning that it is not pos- 
sible to constrain cosmological models using the distance-redshift 
relation in this range. We introduce a simple void parameterisation 
which can reproduce the ACDM distance modulus to sub-percent 
accuracy, and which has a continuous parameterisation from a 
steep and spiky void mimicking ACDM to a smooth one. Thus 
we demonstrate stricto sensu that we can't differentiate between 



A CDM and LTB voids using SNIa alone. However, the general gist 
of lClifton et alj ( l2008h is that the void must be very steep to mimic 
ACDM distances closely, and that this is unnatural. By fitting our 
parameterisation to the new data directly we can see if a steep void 
is preferred over a gentle Gaussian profile. If a steep void is pre- 
ferred by the data, then, this would lend weight to ACDM. How- 
ever, to fully break the degeneracy between ACDM and these void 
models another observable is required. 

As suggested by the curvature test above, a good choice 
is H (z). We use the a ges of passively evolving galaxies to do 
this jSimon et alj|200^ . This probe of the expansion rate of the 
Universe is still relatively new, and thus the data currently says 
little of significance. However, it has the benefit of being a rel- 
atively model independent method to reconstruct H{z). In con- 
trast to other works on this subject we don't use tests like the 
BAO and the CMB. The reasoning for this is mainly that they 
are perturbative tests of the models and so test how perturbations 
evolve; this is particularly important for the BAO and small-Z CMB. 
The theory for this has not yet been worked out so we can't say 
whether results we would obtain have any meaning. Although the 
'background' part of the effect in LT B may be taken into account 
using th e two Hubble rates in LTB jGarcia-Bellido & Haugbjj 



jGarcia-Belli do & Haugb0llel l2008bl ; IZibin et alj I2OO8I: 



Boleiko & W vithe 200^ it has been assumed in previous works 
that the perturbations evolve as they do in FLRW - without any 
scale dependence in their late-time evolution - and so a comov- 
ing sphere at last scattering remains so at late times, modulo the 
distortion from the different Hubble rates in the radial and angu- 
lar directions. However, background curvature enters the Bardeen 
equation for the gravitational potential, and this can vary signifi- 
cantly over a sphere of radius 150 Mpc in void models; this will 
add an important additional distortion which is not yet known, and 
may affect the BAO significantly. 

F urthermore, in LTB models, per turbations are compli- 
cated l lClarkson, Clifton & Febniarvll2009l) because density pertur- 
bations couple to vector and tensor d.o.f.; how important this is in 
dissipating density fluctuations is also not known. Finally, we don't 
consider CMB constraints partly for the reasons discussed for the 
BAO: the large-scale CMB has not been calculated. The small- 
scale CMB may be estimated, however, and intriguingly seems 
to fa vour a non-zero b ang-time function and asymptotic curva- 
ture dClifton et al]|2009h . Though it is not yet clear how degener- 
ate this result might be with the primordial power spectrum, which 
might be important because we don't have an inflationary model 
for void models, this is a very interesting result. Here, however, we 
are really mainly concerned with what the local data tells us about 
the shape of the void, assuming that these voids can evolve out of 
perturbed FLRW, and we don't consider CMB constraints further. 

The paper is organized as follows. In ij2] we describe the var- 
ious void models based on the Lemaitre-Tolman-Bondi metric and 
the kinematical quantities associated with them. Then in ij3] we 
discuss the cosmological data that we have used to constrain the 
models presented in ij2] !j4]is devoted to the data analysis itself, 
and the interpretation of the results; in particular, we discuss the 
ability of data to constrain the smoothness and the size of the void, 
and a series of non-concordance diagnoses. Finally, we conclude 
in ij5]with a summary of the main results and a sketch of future 
developments. 
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2 VOIDS 

2.1 Lemaitre-Tolman-Bondi models 

We model the observable universe as an inhomogeneous void cen- 
tered around us via the spherically symmetric LTB model with met- 
ric 

2 2 ^11 2 2 2 2 

ds — —dt H — , . „ dr +aj_(t,r)r dQ , (2) 

1 — K{r)r'' 

where the radial (a||) and angular (a±) scale factors are related by 

ay = (a±r)' (3) 

and a prime denotes partial derivative with respect to coordinate 
distance r. The curvature k = Hi{r) is not constant but is instead 
a free function. We choose coordinates such that the angular scale 
factor is constant and a±{to,r) — 1. From these two scale factors 
we define two Hubble rates: 



Hi 



a± 

a I 



H\\{t,r) 



(4) 



where an over-dot denotes partial differentiation with respect to t. 
We denote their values today by H±q = _ffx„ (r) = H± {to, r) etc. 
The analogue of the Friedmann equation in this space-time is then 
given by 



„2 _ M K 

^ ~3 2"' 

a I a , 



(5) 



where M — M{r) is another free function of r, and the locally 
measured energy density is 



87rGp(i,r)- 2 

which obeys the conservation equation 
p-f (2ffx + H||)p = 0. 



(6) 



(7) 



The acceleration equations in the transverse and radial directions 
are 



ox 



M iiii 
and 



(8) 



We introduce dimensionless density parameters for the CDM and 
curvature, by analogy with the FLRW models: 



M 



//2 ' 

-l-n 



(9) 



using which, the Friedmann equation takes on its familiar form: 



HI 



(10) 



so Qrn{r) + ^k{r) = 1. Integrating the Friedmann equation from 
the time of the big bang Ib = ts(r) to some later time t yields the 
age of the universe at a given (t, r): 



T{t,r) = t — tB = 



H±,Xr) Jo 



dx 



(11) 



(This integral can be given in closed form, but it's rather ridicu- 
lous.) We now have two free functional d.o.f.: Q.m{r) and tB{r), 
which can be specified as we wish. However, if the bang time 
funct i on is not c onstant this represents a decaying mode jSilkl 
1 19771 : 1 Zibinll20o3) : consequently we set ts = throughout, which 



means that our model evolves from FLRW. As a result, the age of 
the universe r is then a constant, and equal to the time today to. 
So, by solving Jilt for H±„ (r), we have that: 



2 

— Qm sinh"-^ J ^ 



Hfc < 

fife = (12) 
fifc > 



When we introduce the Hubble constant below, we shall use Ho = 
H±g(r = 0), which fixes to in terms of Ho, ^m{r = 0) and 
^kir = 0). 



2.2 Distance modulus and observables 

In LTB models, there are several approaches to finding obs ervables 
such a s distances as a function of redshift. We refer to lEngvislI 
( l2008h for details of the approach we use here. On the past light 
cone a central observer may write the t, r coordinates as functions 
of z. These functions are determined by the system of differential 
equations 



dt 

dz 

dr 
dz 



{l + z)a\\H\\' 



(13) 



(14) 



where (t, r) = H\\{t[z),r{z)) — _ff|| (2), etc. The area distance 
is given by 



dAiz) = a±{t{z),r{z))r{z) 



(15) 



and the luminosity distance is, as usual dL{z) = (1 + z^dA^z). 
Other observables can be calculated from these relations; in partic- 
ular, the distance modulus is given by 



^{z) — m — M — 5 logj^ 



dL(z) 



IMpc 



+ 25, 



(16) 



where m is the apparent magnitude of a source with intrinsic mag- 
nitude M. The relative ages of galaxies may be calculated from 
Eq. ([13} to give /fy (z). 

The algorithm for calculating functions of z is as follows: 

(i) Choose a profile Sim. (?•) today; this gives ^Ik (r) and H±g (r). 

(ii) At any spacetime point {t, r) we have a_i_{t, r) by inverting 
Eq. (TTJ- 

(iii) We may then calculate ay {t, r), and hence iify (t, r), from 
Eq.©. 

(iv) Integrate Eqs. il3l and J14t from z = with initial con- 
ditions t{z = 0) = to and r{z = 0) = to give the parametric 
equations t{z),r{z). 

(v) Calculate all functions as a function of z. 



2.3 Void Profiles 

For our main model we introduce a profile which is capable of re- 
producing the ACDM distance modulus to high accuracy, as well 
as being able to control the smoothness at the centre of the void via 
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Figure 1. Void #1 shown for different values of v. For ^ — > oo we liave 
a sharp, non-differential void which can mimic the ACDM distance mod- 
ulus at low z (middle right). (The effective deceleration parameter, shown 
bottom, is defined later) Only when the profile has discontinuous derivative 
(middle left) can the deceleration parameter be negative at 2 = 0. 



the parameter v. The parameterization is given by (for u > Q and 

r^out ^ir 



V tanh ■ 



' tanh — , 
a J 

(Model # 1) 



(17) 



where f2i„ and fiout are the value of Q,m at the centre of the void at 
infinity, respectively. The parameter a characterizes the size of the 
void, but a more useful quantity is the full width at half-maximum 
(FWHM), calculated by solving 



for r„ 



1 



(^out H~ ^in 



(18) 



, p^yj^„ numerically. (We assume 0,,^ [r) is even about r = so 
thatn„(-T-) = VL^{r).) 

The parameter v is chosen so that when v is finite this function 
IS at the origin (its first and second derivatives are well-defined 
and equal to zero at r = 0). In the limit — >■ cxd we have 



T 

fim(r) — 5> fiin + (r^out - fiin) tanh — , 



(19) 



which can give an extremely good fit to the ACDM distance mod- 
ulus - at the expense of not being differentiable at the origin. The 
parameter u gives us the power to control the sharpness of the void 
at the origin: the larger v is, the steeper the void; the smaller u is, 
the flatter the void is at the centre. Finally, note that the indeter- 
minate form of Q.m{r) at = 1 is only due to the form of the 
parameterization: we can analytically continue the function (TtJi to 
include v = 1: 



fiin) tanh ■ 



— sech^ 



-1 



(20) 



In Fig.[I] we show this void for different values of i^. 

In addition to the parameterization given by Eq. (TTJ, we also 




Figure 2. For a given flat ACDM model we find the ability of each void 
model to reproduce the same Here we show the largest differences as 
a function of the matter content of the FLRW model for the best-fitting void 
model in each case. 



considered the following void profiles: 

# 2 : Q.m{r) = riout — (fioiit — fiiii) exp 



#3: ^^{r) 
#4: n„(r-) 
#5: n„(r-) 



^OLit (^out ^in) " 



- (a. 

- (a. 



(7) 



Typically we will fix Sloui ~ 1, so that the spacetime is 
asymptotically flat, in keeping with generic predictions from infla- 
tion. However, early-universe models which might produce a void 
of the kind we are co nsidering have not been explored (although 
see lLindeetalJ ( ll995l) ). On the other hand, because we have set the 
bang time to ts = 0, the models we consider evolve from a per- 
turbed FLRW model at early times (the void has \Sp/p\ ~ 10^'^ at 
last scattering), so this may conceivably be natural. 

It is known that there exist LTB models which can 
give t he FLRW distance modulus exactly for a central ob - 
server l lMustapha et al.ll997l : lcilerier et alj2009l : lYoo et al.l2008h . 
Profile #1 is a void parameterisation which can accurately mimic 
ACDM to high precision. If we assume Qout ~ 1 and perform a 
least-squares fit of the void ij,{z) to ACDM models for < 2: < 
1.6, then we find that all of the profiles given can reproduce the 
distance modulus of ACDM to sub-percent accuracy - see Fig. [2] 
Model #1 can produce a ACDM distance modulus to very high ac- 
curacy, which requires v — 00 {a spiky void). The corresponding 
radial profiles for our 5 different best-fitting to ACDM void models, 
as well as their distance moduli and effective deceleration parame- 
ter (defined below by Eq. i24i ) are shown in Fig. [3] For Model #1, 
even though the distance modulus is effectively the same as that 
of FLRW, the deceleration parameter is noticeably different; this is 
because the Hubble rates are different. 

Finally, note that with Sloui = 1, voids 2-5 have 3 parameters 
(including Ho), which is the same as a curved ACDM model, while 
#1 has 4. Note also that the parameter a has dimensions of length, 
an issue we will return to. 
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Figure 3. Present day radial profiles (top 6 panels) and corresponding dif- 
ference in distance moduli for each void and that of ACDM (third row from 
top), as well as the resulting effective deceleration parameter as a function 
of redshift (bottom), after fitting each void with Qout = 1 to ACDM with 
= 0.3. 



2.4 Physical length scales and the distance modulus 

In LTB models, subtleties arise concerning the relation between J/q 
and the magnitude-redshift relation. In the FLRW case, since Ho is 
just a magnitude offset related to the intrinsic (absolute) luminosity 
of a SNIa, it is usually removed from the analysis by marginalizing 
over it. The intrinsic absolute magnitude of SNIa's is poorly con- 
strained, and since this value and the value of are degenerate in 
ACDM, the Hubble constant is poorly measured by SNIa. 

In an LTB model we have two length scales which are inde- 
pendent: the Hubble length, H^^, associated with the expansion 
time, and the void scale depending on the physical 'size' of the 
underdensity. In the models considered here, this may be character- 
ized by the parameter a or the FWHM. The ratio of these two scales 
enters as a dimensionless number, and this must therefore have a 
physical significance. Thus, when we change Ho keeping the void 
scale fixed this will be reflected in the distance modulus in a non- 
trivial way, namely, if one rescales the distance by r i— >■ r/Ho, in 
the void profiles this becomes r/a = r/aHo. 

As a result, higher (lower) values of Ho not only shifts the 



Figure 4. Plots of - ^^p'S" for FLRW and LTB models for different 
choices of Hq- This illustrates that, in the LTB case, Hq not only affects 
the vertical displacement of fJ,{z), but also its shape. This means that our 
obtained void parameters are partially dependent on this value, and thus a 
standard normalisation cannot be applied. 



distance modulus curves down (up), but it also affects the overall 
shape too. This is shown in Fig. [4] where we can see that the shape 
of the fj,{z) function changes with Ho- If crHo were held fixed when 
Ho is changed (for example by fixing a in units of h^^Mpc rather 
than Mpc, say), this would compensate for this effect and then Ho 
would be a pure normalisation of the distances. However, linking 
the two independent scales in this way is rather restrictive. 

This is an important issue, and, as far as we are aware, has 
not been previously considered. This means that the shape and size 
of our best fit void model is dependent on the value of Ho we ob- 
tain when we fit the voids to data. With SNIa being able to poorly 
constrain the value of Ho, this leads to an additional uncertainty 
in the best-fitting parameters obtained. The best-fitting value of Ho 
as indicated by the latest supernova dataset is clearly not in agree- 
ment with other measurements, such as that found by the HST Key 
Project using Cepheid variables. However, we estimate that the ad- 
ditional uncertainty in void models is below the 5% level, and as 
such does not play a significant role in the current error budget. Fu- 
ture studies of inhomogeneous models, be it globally LTB models 
as in our case, or mass-compensated ones in other cases (see e.g. 
iKrasihski & Hellab\] ( l2004h . [Alexander et alj | |2007|) ). wifl need to 
consider the effect of Ho on their results, possibly by fitting the 
supernova light-curves simultaneously with the model. 



3 COSMOLOGICAL DATA 

In this paper, we confront each of the 5 void models we introduced 
in the last section with the largest sample of SNIa to date; the C on- 
stitution dataset consist ing of 397 SNIa Jllicken et al.ll2009bl) . as 
well as the H{z) data l lSimon et al.ll2005l) . consisting of 10 points 
when we include the HST Key Pr oject value of the Hub ble con- 
stant, //o = 72 ± 8 kms^^Mpc"^ jFreedman et alj200lh. 

The Constitution dataset of SNIa jHicken et aU l2009bl) 
is com prised of a l arge s ample of nearby (z < 0.08) ob- 
iects |Hicken et aU l2009al) combined with the Union dataset 



jKowalski et al.l 20081) which covers a redshift range up to z 
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1.85. This sample comprises 397 spectroscopically confirmed 
SNIa, making it the largest publicly available and uniformly anal- 
ysed SNIa dataset to date. An additional intrinsic dispersion of 
= 0. 12 is added to the errors of each supernova to better estimate the 
error due to the uncertainty in absolute magnitude of these events. 

In ^ we consider several other supernova datasets. These 
samples do not include the large number of SNIa's at low redshift 
that comprise a lar ge proportion of the Constitution dataset. The 
Davis et al. dataset jPavis et aLl2007l) combines SNIa's discovered 
from the ESSENCE (Equation of State: SupErNovae trace Cosmic 
Expansion) survey with those from the Hig h-Z release to produce 
a sample of 192 objects. The Union dataset l lKowalski et"al]l2008l) 
comprises 307 objects that have been combined fro m other pub licly 
released datasets, whilst the ConstitutionT sample l lWeil2009l) uses 
the Constitution dataset described above, with 34 SNIa removed in 
order to make the sample be in better agreement with other SNIa 
datasets (not to mention the ACDM model!). 

The expansion rate of the Universe, H{z) has been con- 
strained using the differential ages of passively evolving galaxies 
as determined by fitting SED te mplates to their spect ra. Using data 
from the Gemini Deep Survey ( Abraham et al. '2004") and a sample 
of field early-type galaxi es jTreu et alJll999. 2001, 20021) along 
with two radio ga l axies ([Dunlop et al.lll996l : ISpinrad et al.|[l997t 
iNolan et al. 

Illooi), I 

Simon et al. I( l2005l) were able to constrain the 
evolution of H{z) for 0.1 < z < 1.8. In this analysis we in- 
clu de the measurement of the local value of H{z) as determined 
by iFreedman et al.l 1 I2OOII) to make 10 data points, but note t hat a 
recent, more accurate m easurement of Hp jRiess et aL 1 1200^ and 
additional H{z) points jStem et al]|2009h have since been deter- 
mined, and would place somewhat tighter contraints on the models 
considered in this work. 

Thus this analysis uses two independent probes of the expan - 
sion history of the Universe. Other work, such as IZibinI ( 1200 Sl) . 
combine supernova datasets with other probes, such as BAOs and 
CMB measurements; however, these are not considered in this anal- 
ysis, as explained in 21 since the growth of perturbations has not 
been properly explored in UTB cosmologies and thus approxima- 
tions to ACDM relations may not be valid. 



4 DATA ANALYSIS 

4.1 Overview of Numerical and Statistical Method 

We made us e of the publically available easyUTB Fortran90 code 
provided by iGarcia-Bellido & Haugb0li3 ( l2008ah to run through 
more than a million different parameter values, simultaneously 
computing iJ.{z), H{z) and the statistic for each model with the 
intent of finding the best-fitting parameters for the voids consid- 
ered. The code is setup to compute for the SNIa (Xsjv/a)' BAO 
(Xb_4o) and CMB data (Xcj/b)' ^^''^ '^^^ testing of the 

UTB models against the last two for future work due to the fact that 
these quantities are affected by the growth of structure in the uni- 
verse. On the other hand, we have added the Hubble rate data into 
our analysis, and we have us ed the longitudinal Hubble rate i^y as 
our model prediction, since jSimon et al.ll200^ determined H{z) 
via dt/dz measurements, which in the UTB case (see Eq. l |13t ) is 
related to H'^\. The main reason for adding the H{z) data into our 
analysis is because, as shown in ij2l all of our void parameteriza- 
tions can mimic /i^*^"'", making it impossible for us to distinguish 
between two such models if the data prefers ACDM; so using H{z) 
in addition to the SNIa data (which we refer to as SNIa+H) allows 




r(Gpcl 
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Figure 5. Best-fitting Q^ii") profiles for void model #1 from various SNIa 
datasets. 



us the further opportunity to find detectable differences between 
ACDM and UTB void models. The easyUTB code is not set up to 
compute models in which ilrn{r) > 1. For this analysis, we mod- 
ified the code to allow ilm{r) to take on any value, as required by 
our oscillating voids, although for our main investigation we have 
fixed Slout to 1 . 

Regarding the fitting of our void models to data, we computed 
the x^ statistic, and associated reduced x'^, Xred' to determine the 
goodness of fit of the model to the data. Cases where Xrcd < 1 
indicate that the additional intrinsic dispersion added to the type 
la supernova error estimates is a conservative choice. On top of 
this, the confidence limits on each of the parameters was calculated 
using their likelihood distributions. 



4.2 Pointy or Smooth? 

As demonstrated earlier, void models that are pointy at the origin 
are capable of reproducing fi^^^'^ to arbitrary accuracy. However, 
what level of smoothness/sharpness at the origin do the SNIa data 
demand? 

In Table[T] we show the best-fitting model parameters of void 
#1 whe n fitting to the sam ples described in ^ T he Davis et al. 
samp le jPavis et al.l2007l) (A), the Union sa mple 1 Kowalski et al.l 
l2008h (B), the original Constitution sample jHicken et al]|2009bh 
(C), the Constitution sample with intrinsic dispersion added (D , our 
princ iple SNIa dataset) and finally the ConstitutionT sample l lWeil 
llooi) (E) are considered. From these results we see that the best- 
fitting u is different for different data sets. The reason is as fol- 
lows: in the older data sets such as A and B, there are few to none 
low-redshift SNIa, so that f is essentially meaningless, although, as 
shown in Table[T]and Fig. [5] the data sets given by A and B choose 
sharp voids naturally. On the other hand, with the abundant low- 
redshift SNIa in data sets C-E, 1/ can play an active role in determin- 
ing the gradient of the profile around the origin, and, in particular, 
these data sets favour u ~ 1, corresponding to a voi d that is very 
smooth around the origin. This supports the claim bv lClifton et al.l 
120(3) that "realistic" voids should be smooth around the origin. 
Moreover, as Fig. |6] shows, void model #1 turns out to be flatter 
around the origin than the Gaussian model (#2). Note that even with 
the next generation of SNIa measurements, a very low z distinction 
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Table 1. Best-fitting parameters of void model #1 from various SNIa data sets. 



Data Ho Un, a log^QU XsNia xL FWHM 
Set # Gpc Gpc 



A 


65.78 


0.065 


9.26 


2.95 


195.38 


1.03 


10.19 


B 


70.60 


0.095 


6.50 


5.00 


310.68 


1.02 


7.14 


C 


64.18 


0.140 


2.17 


8xl0-'3 


459.93 


1.17 


5.44 


D 


64.17 


0.149 


2.05 


4x10"'' 


308.82 


0.78 


5.13 


E 


64.00 


0.158 


2.33 


3.4x10-2 


270.88 


0.75 


5.52 



£2 (r) 



Hir) 




k(r)/lk(0)l 








-0.2 




-0.4 




#5 


-0.6 




#4 






#3 


-0.8 




#2 


-1 




#1 



Zooming ill around ihe origin 

-0.7 \ 



5 

r(Gpc) 




Figure 6. Best-fitting void profiles from SNIa data. 
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Figure 7. Best-fitting void profiles from SNIa+H data. 



Table 3. Best-fitting parameters for the flat ACDM model from SNIa and 
SNIa-l-H data. See the caveats in i|2.4l regarding what the value of Hq means. 



Data 


Ho 






■^red 


SNIa 


64.93 


0.29 


311.18 


0.786 


SNIa+H 


64.55 


0.32 


326.51 


0.804 



will not be realistic, since below z ~ 0.015 SNIa are dominated 
by peculiar velocities and are not properly in the Hubble flow - 
indeed, a cosmological model becomes a bit meaningless on such 
small scales. 

In Table|2]we show the best-fitting parameters for void models 
#2-5 after fitting them to the Constitution SNIa sample with intrin- 
sic dispersion added (D). It is interesting to note that all of these 
models fit the SNIa data equally well compared to void model #1 - 
see Table [T] In addition, the FWHM of each void is roughly the 
same. Therefore, the simple void profiles such as the ones we have 
studied can all agree on the characteristic size that a void should be 
to fit the data, implying that one could effectively use any of these 
toy models in order provide a reasonable estimate of the physical 
size a typical void should be in order to fit the data. Comparing 
the x^'s of the voids (model #1 included) with that of the ACDM 
model in Table[3]we obtain a similar (although slightly better) fit to 
the SNIa data as flat ACDM. 

Given the difficulty of distinguishing LTB voids from ACDM 
using the SNIa data alone, we performed a combined analysis of 
the SNIa sample denoted D above with the H{z) data. In Table |4] 



we show the best-fitting parameters from SNIa-l-H data for each of 
our 5 void models, along with the corresponding x^'s and FWHM. 
Notice again that the x^, and Xred values for the voids are not only 
comparable to that of ACDM (see Table O, but are slightly lower. 
In Fig.|7]we show the resulting void profiles for the best-fitting to 
SNIa-l-H data. The sizes (FWHM) of the voids are slightly larger 
(by roughly 1 Gpc) in this case compared to that of the fit to SNIa 
case. This is a result of the fact that the fit to H{z) by itself favours 
enormous voids (FWHM ~ 65 Gpc!), and thus when combined to 
the SNIa data, we obtain bigger voids. 

In Fig. [8] we show the best-fitting to SNIa-l-H data residuals 
(top panel) and Hubble rate (bottom panel) for each of our 5 differ- 
ent void models, as well as that of ACDM and EdS. In each case, 
the respective data points are overplotted (grey circles), along with 
I-cr error bars. To illustrate the difficulty in distinguishing these 
voids from that of ACDM using SNIa data alone, in Fig.|9]we plot 
the magnitude difference between the two models after fitting our 
5 different voids to the SNIa-l-H data. 

4.3 Likelihood Contours 

In this section, we explore the degeneracies between all the possi- 
ble pairs of parameters by constructing joint-parameter likelihood 
contours, and also determine the likelihood distributions of each 
parameter from SNIa and SNIa-FH constraints for void model #1. 

Table |5] shows the marginalized best-fitting parameters for 
void model #1 with 95 per cent confidence limits, from fitting to 
SNIa and SNIa-l-H data, where the priors are also given. We see 
again that the main difference between the parameter values in the 
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Table 2. Best-fitting void model parameters from SNla data. 



Model Ha U,, a x%Nla xL 

# kms^^Mpc^^ Gpc Gpc 





64.36 


0.140 


3.29 


308.83 


0.782 


5.47 


3 


64.39 


0.129 


3.21 


309.36 


0.783 


6.42 


4 


64.34 


0.148 


1.35 


308.54 


0.781 


5.77 


5 


64.35 


0.143 


1.90 


308.70 


0.782 


5.29 



Table 4. Best-fitting void model parameters from SNla+H data. 



Model* Ho n,„ a logio xld FWHM 



1 


64.24 


0.119 


4.94 


0.71 


321.64 


0.796 


6.81 




64.33 


0.120 


3.77 




321.74 


0.794 


6.28 


3 


64.38 


0.114 


3.59 




321.57 


0.794 


7.17 


4 


64.30 


0.125 


1.55 




321.81 


0.795 


6.65 


5 


64.35 


0.121 


2.21 




321.86 


0.795 


6.14 



SNIa and SNIa-l-H cases is the value of a: the SNIa-l-H case yields 
a higher value for a, which arises from the fact that, when fitting 
void models to H{z) data only, extremely large {a ~ 20) Gpc 
and empty (fiin ~ 0.01) voids are the ones that give the minimum 
X^, so when combining with the fit to SNIa results which has a 
lower value for a, it is only natural that a is then larger in the fit to 
SNIa-l-H case than in the fit to SNIa case only. Notice in the column 
for the smoothness parameter v (or logioi^ as shown), that the best- 
fitting values are roughly 3 orders of magnitude larger than that in 
the unmarginalized case (see Tables [2] and |4j, indicative of sharp, 
spikier voids. Also note that there is no upper bound on only a 
lower one, but we shall return to this topic in a moment. 

In Fig. [To] we show the joint-parameter likelihood plots from 
SNIa-l-H constraints. The inner, red-filled regions represent the 68.3 
per cent confidence level, the orange regions 95.4 per cent, and 
the off-yellow regions correspond to the 99.7 per cent level. The 
crosses indicate the positions of: the best-fitting to SNIa data (blue), 
the best-fitting to H{z) data (green), the best-fitting to SNIa-l-H data 
(cyan), and the best-fitting marginalized (black), parameters. 

Let us quickly discuss the degeneracies that exist between the 
parameters. The strongest degeneracies are between all the possible 
pairs of Q.-^, cr and v, and can be interpreted as follows: in order 
to give the same x^, if one wants a larger void (larger cr), then 
must decrease , and the void must be sharper at the origin (larger p), 
and vice versa. The degeneracies between those same parameters 
and Ho are less obvious, but still present: emptier and sharper voids 
require a larger Ho, whereas a larger void needs a lower Ho in order 
to give similar x^, and vice versa. 

In Fig.[TT]we show the one-parameter likelihood distributions 
for each parameter from SNIa-l-H constraints. The vertical lines 
correspond to: the best-fitting to SNIa (blue, space-dashed), best- 
fitting to H{z) (green, dot-dashed), best-fitting to SNIa-l-H (cyan, 
solid) and marginalized best-fitting (black, dashed), parameters. 
Notice that Ho and Slj,, are fairly well constrained, whereas cr and, 
even more so, v, are very poorly constrained individually owing to a 
degeneracy between them; taken together, however, the constraints 
in the cr — plane are reasonable. Note that v is not constrained 
from above because = oo is the ACDM-void by construction. 



Table 6. Information criteiion for the best-fitting void models and that of 
ACDM model from SNIa-l-H constraints 



Model 


AlC 


BIC 


1 


329.70 


345.68 


2 


327.78 


339.77 


3 


327.60 


339.60 


4 


327.84 


339.84 


5 


327.89 


339.89 


ACDM 


330.52 


338.53 



4.4 Which model is most preferred? 

According to the Xred statistic, the void models considered here are, 
very slightly, favoured over ACDM. However, there are more so- 
phisticated techniques applicable to model sel ection, such a s the 
(corrected) Akaike information criterion (AIC) ( lAkaikeiri974l) 

AIC = xL + 2fc + ^^^^, (21) 

where k is the number of parameters and A*" the numbe r of data 
point s, and the Bayesian information criterion (BIC) l lSchwarzl 
[l97i) 

BIC = xL + fclnAr. (22) 

A lower AIC or BIC implies a more favoured model. In Table[6]we 
summarize the AIC and BIC values for each void model as well as 
that of ACDM. Using the AIC criteria we can see that all models 
are more-or-less equally favoured, whilst the BIC, which heavily 
penalises additional parameters, moderately disfavours void model 
#1, in comparison to the remaining void models and ACDM. 



4.5 Non-Concordance Diagnostics 

Finally, we consider quantities that can serve as useful conceptual 
tools for distinguishing FLRW/ACDM models from LTB models. 
Because the voids are really just toy models, constructing diagnos- 
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Table 5. Marginalized best-fitting parameters at 95 per cent confidence from SNIa and SNIa+H constraints. Note 
the amusing error margins on the last two parameters individually ! The error in the u — v projection are much 
more informative - see Fig. (To) 



Data Hq f2i„ a (Gpc) logio ^ 



Priors 


58.00- 


74.00 


0.00-0.24 


0.00-^ 


-30.00 


-1.30 


-3.70 


SNIa 


64.61 


f 1.69 
-1.39 




6.72 


f 17.68 
-5.47 


3.3l 


OO 

3.97 


SNla+H 


64.72 


f 1.58 
-1.49 




7.84 


f 12.54 
-6.27 


3.5l 


oc 
3.84 



Distance Modulus 



Distance Modulus 
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Figure 8. Top panel: Best-fitting SNla-HH distance moduli, normalized to 
that of a Milne universe with Ho given by that of the best-fitting ACDM 
model, for each void, along with that of the best-fitting ACDM and EdS 
models, as well as the Constitution SNIa data points (grey circles). Bottom 
panel: The best-fitting SNIa+H Hubble rates for each void model and that 
of ACDM and EdS, with the H{z) data points overplotted (grey circles). 
Notice how well the void H{z) curves follow the data points. 



tics from ACDM models allow us to visualize what our real con- 
straints on ACDM are. 

Consider the standard deceleration parameter in FLRW mod- 
els which depends purely on the Hubble rate and the observable 



Figure 9. Difference in distance moduli between the best-fitting SNIa-HH 
voids and that of ACDM. Up till a redshift 2: ~ 1, these void models are 
practically indistinguishable from that of ACDM. There is no way to tell 
the difference between ACDM and our void models from SNIa alone, but 
as the bottom panel of Fig.[8]shows, H{z) may play a key role in the future 
for testing for deviations from ACDM. 



redshift: 



q{z) 



-1 + 



(1 + ^) d 
H{z) dz 



H{z) 



(23) 



In testing for acceleration in FLRW models, we would like to ide- 
ally c onstrain this function directly from data jMortsell & Clarksonl 
l2009h . Now, using Eq. fTsl l. we can define an effective deceleration 
parameter in void models as 



Cff / N 
1 (^) 



-1 



{l + z) d 



Hiiiz) 



(24) 



dz' 

We anticipate, given the results of I Vanderveld et al] ( |2006|) . that 
models containing a weak-singularity would have qf)^ < 0, and 
models without the singularity would give q^^ > 0. In Fig.ll2lwe 
show the effective q'^^^z). These all have qo > w hich is expected 
for smooth voids. It is interesting to note that IShafieloo et al.l 
( I2OO9) find go > for various classes of dark energy models within 
the FLRW paradigm when using the same datasets we consider 
here. In fact, it is worth comparing our q''^\z) with the form they 
find when they fit only SNIa - they are qualitatively very similar. 

We may also consider an effective dark energy equation of 
state for the void models which we define as 

2(1 + z)d'J + 3< 



eff / N 



3[H§n^{l + Z)^d',^ -l]d'. 



(25) 
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Figure 10. Joint-parameter likelihood contours from SNIa+H constraints. 
Black crosses indicates the marginalized best-fitting values, blue crosses 
show the best-fitting to SNIa values, green crosses are the best-fitting to 
H(z) values, and cyan crosses shows the best-fitting to SNIa-l-H values. 
See ^2.4l for a discussion of what fits to Ho mean. 



where dc = (1 + z)dA is the comoving angular diameter distance. 
This is the EOS which, in a flat FLRW dark energy model with 
energy density Q.m gives the comoving angular diameter distance 
dc(z)- This gives us another nice way to visuahze the differences 
between our voids and more standard dark energy functions used 
in the standard model. We show in Fig. [13] this function for our 
best-fitting void models. Note the apparent phantom behaviour for 
moderate redshift. 

Finally we consider two tests which have recently been pro- 
posed. The first is a test for the flat ACDM scenario: 



(26) 



should be constant if the concordance model is cor 



rect 



Zunckel & ClarksonI |2008|) (called Om{z) in ISahni et al.l 



( l2008|y K which implies that 

^{z) = 2[{l + zf -l]Hod'J{z) 

+3(1 + zfHod[{z) [1 - {Hod[{z)f] , (27) 

is zero in the concordance model jZunckel & Clarksoij|2008l) . 

The second test is a much more general test of a FLRW geom- 




Figure 11. ID likelihood distributions for each parameter from SNIa+H 
constraints. Black dashed lines corresponds to marginalized best-fitting val- 
ues, while the blue, green and cyan lines con'espond to the best-fitting to 
SNIa, H(z) and SNIa+H values, respectively. 
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Figure 12. Effective deceleration parameter for the best-fitting models to 
SNIa+H data for our 5 different void models, as well as that for ACDM and 
EdS. Note that for all these voids models we have an effective deceleration 
today. 



etrv dClarkson et alj200§) . 



fife yz) = 



ld,{z)HQf 

should be constant for any FLRW dark energy model, and 

^(^) = l + H^^{zf[d,{z)d['{z)^dUzf] 
+H\\{z)H\\{z)d,{z)d[{z) , 



(28) 



(29) 



should be zero. The utility of these tests is that they can be used 
independently of any model at all. Considering these functions in 
our void models helps us visualize the difference from FLRW in 
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Figure 13. Effective equation of state parameter from distance measure- 
ments for tlie best-fitting models to SNIa-HH data for our 5 different void 
models, as well as that for ACDM and EdS. 
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Figure 14. Top: Effective n^^iz) and ^^=^,''(2) for our 5 best-fitting void 
models to SNIa-l-H data, along with that of the ACDM model. For the case 
of 0^(2), we have placed the H{z] data used in fitting the models, without 
the error bars (which are rather large). Bottom: the corresponding second- 
derivative tests Jf{z) and ^{z). 



terms of observable functions. This is shown in Fig. [14] Note that 
for Q,'^{z), our best-fit ting voids produce a curve which is very 
similar to that found bv lShafieloo et al.l ( |2009|) considering FLRW 
dark energy models. 



5 DISCUSSION 

In this work, we have discussed LTB void models as solutions to the 
dark energy problem, and have considered how to constrain them 
with cosmological data. Of particular interest is the steepness of the 



void near the centre of symmetry, as well as the need to combine 
independent data sets in order to distinguish voids from ACDM. 

As far as SNIa and distance data are concerned, it is often 
claimed that to fit the data a spiky void near the centre is required; 
this is not the case. To reproduce the ACDM distance modulus ex- 
actly, one needs a void with a point at the origin; if one fits the ac- 
tual data, and not a prediction of the ACDM model, smooth voids 
are perfectly compatible with observations, and their likelihoods 
make them indistinguishable from a ACDM scenario. By introduc- 
ing a new parameterization we can smoothly parameterize away 
from the "ACDM-voids", to test specifically how steep the void 
would have to be. 

Moreover, we have included H{z) data in the hope of remov- 
ing the degeneracy between a class of voids and the ACDM mod- 
els. Unfortunately, current uncertainties on the age measurements 
do not yet allow for such a discrimination between the models, with 
only the void model with an additional parameter being moderately 
disfavoured by the BIC compared to ACDM. All other models are 
equally favoured by current datasets. 

On the other hand, a careful look at the joint constraints from 
SNIa and H{z) data provides evidence for a robust estimation of 
the void size through the full width at half maximum. Indeed, in 
agreement with previous results, we have found that the size of 
the void is of order 6 Gpc. This estimation appears to be fairly 
independent on the parameterization of the void profile. 

Finally, we have illustrated how our void models look when 
viewed in terms of standard FLRW functions. For example, the 
voids, interpreted through the lens of a flat dark energy fluid give a 
'w{z) which has phantom behaviour at intermediate redshifts, and 
a deceleration parameter which is positive for low redshift. These 
profiles may be usefully considered as signatures of voids. Further- 
more, we have also shown profiles for iT^(z), which is constant 
only for flat ACDM, and flf{z) which is constant in any FLRW 
model whatsoever. These will be very useful diagnostics in future 
studies of dark energy because they allow us to measure deviations 
from flat ACDM and FLRW models respectively in a model in- 
dependent way. By seeing what they look like in void models, we 
understand further the range of non-concordance behaviour possi- 
ble in these diagnostics. It is rather interesting to note, in fact, that 
the effective deceleration parameter, as well as the effective mat- 
ter density predicted by LTB models look very similar to previous 
estimations of these quantites rel ying on a parameteriz ation of the 
equation of state of Dark Energy jShafieloo et alj200^ . 

As a final remark, re cent constraints on no n-standard cos- 
mologies were presented in lSollerman et al. using the first- 

i ear Sloan Digital Sky Survey-II (SDSS-II) supernova results 
Frieman et alj|2008n in combination with other data, such as the 
latest BAO and CMB measurements. In p articular, using a simple 
Gaussian void profile allowing fiout to varv. lSollerman et al] l l2009h 
demonstrate that with this new dataset LTB models are still capable 
of fitting the S NIa-only data be t ter tha n the ACDM model. 

However. lSollerman et al.l ( [2009h do not find evidence that the 
additional d.o.f. they allow when compared to our Gaussian profile 
is supported by the latest dataset. This dataset contains new SNIa, 
which now fill in the gap around z ~ 0.2, but lack the recent low 
redshift supemovae that are included in the Constitution dataset. 
We note that their conclusion that LTB models are a worse fit than 
ACDM are as a result of using BAO/CMB data which we do not 
use here for the reasons discussed above. Thus, we expect that their 
new SNIa data will affect our conclusions only marginally. 

This work strongly suggests the importance of testing LTB 
models with further observations such as CMB and BAO, jointly 
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with SNIa and age measurements, in order to discriminate between 
models. Such analyses, however, will rely on a careful implemen- 
tation of structure for mation and CMB emission in an LTB back- 
ground as initiated in IClarkson. Clifton & Februarvl ( l2009h . This 
will be the topic of forthcoming studies. 
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